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ABSTRACT 

Giant Lya nebulae (or Lya "blobs") are likely sites of ongoing massive galaxy formation, but the 
rarity of these powerful sources has made it difficult to form a coherent picture of their properties, 
ionization mechanisms, and space density. Systematic narrow-band Lya nebula surveys are ongoing, 
but the small redshift range covered and the observational expense limit the comoving volume that 
can be probed by even the largest of these surveys and pose a significant problem when searching 
for such rare sources. We have developed a systematic search technique designed to find large Lya 
nebulae at 2 < z < 3 within deep broad-band imaging and have carried out a survey of the 9.4 square 
degree NOAO Deep Wide-Field Survey (NDWFS) Bootes field. With a total survey comoving volume 
of ^10^ h^Q^ Mpc^, this is the largest volume survey for Lya nebulae ever undertaken. In this first 
paper in the series, we present the details of the survey design and a systematically-selected sample 
of 79 candidates, which includes one previously discovered Lya nebula. 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: high-redshift — galaxies: 
surveys 



1. INTRODUCTION 

Gia nt radio-quiet Ly a ne bul ae fe.ps;.. iFrancis et al.l 
19961: iS teidel et al.' "2000" 'Palunas et al." "2004; 
Matsuda et al. 2004; Dev et al. 2005; Saito et al. 2006\ 
Nilsson et al. 2006; Smit h Jarvisi i2007( : iGreve et al. 
2007; Yang et al. 2009) likely provide an observational 
window into the physics of ongoing massive galaxy 
formation. As large as ^100 kpc across, Lya nebulae 
emit copious Lya emission (~ 10^^ erg s~^), signaling 
the presence of highly energetic phenomena. Lya 
nebulae are sometimes surrounded by or associated with 
young, star -forming galaxy popu l ations and obscured 
AGN (e.g., 'Cha pman etldl [200TI: iBasu-Zych &: Sd^ 
[2004; Matsuda et al. 200l; iDev et al.l [20051 : iGeach et al.l 
[2007> . .2009l Prescott et al. 2012c, submitted) and 
have been show n to res ide in over d ense envir onments 
(iMatsud a erall [2004 120051 "2009; Sait o et al.l 120061 : 
[Prescott ~al.l 120081 : lYang et a l. 2009, 20 ij). Studying 
Lya nebulae can provide important insight into the 
process of galaxy formation, but after a decade of study, 
we have not yet developed a complete understanding of 
the dominant power source or, even more fundamental, 
the space density of Lya nebulae. 

The problem lies in having such a small sample - 
only a few dozen large 100 kpc) Lya nebulae are 
currently known - and in attempting to measure the 
space density of a rare population over restricted vol- 
umes. Radio-quiet Lya nebulae were first discovered in 
targeted narrow-ban d imagin g studies of known over- 
densities (e.^.. Fra ncis et al.l J, 996; Stei del et al.l [20001 : 
IMatsuda et al.ll2004f ). Since then, the sample has slowly 
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grown as a result of ever more expensive narrow-band 
surveys over increasingly wide areas (e.g., Nilsson e t al. 
[2006 :^ Smith fc Jarvis 2007: Saito et al. 2006; Yang et al. 
[20091 . [20101 : [Matsuda et atllloTl . Building a sufficient 
sample of such rare objects over a wide range of envi- 
ronments and accurately measuring their space density 
requires surveying larger comoving volumes than are typ- 
ically feasible via standard narrow-band imaging. 

We have therefore taken a complementary approach 
by developing a systematic search algorithm to find lu- 
minous Lya nebulae at 2 < z < 3 using deep broad-band 
imaging and have tested the technique using archival 
imaging d ata from the NOAO D eep Wide-Field Survey 
(NDWFS; iJannuzi fc Dev[[T999[ ). We describe our sys- 
tematic broad-band Lya nebula survey in three papers. 
This paper (Paper I) discusses the survey design and 
candidate selection approach. In Section [21 we describe 
the algorithm used to select Lya nebulae. Section [3] 
presents the sample of 39 first priority and 40 second 
priority Lya nebula candidates selected from archival 
imaging of the Bootes Field of NDWFS and discusses 
the potential sources of contamination. In Section HI 
we address the complementarity of our approach relative 
to traditional narrow-band surveys, and we conclude in 
Section m Discussions of (a) our successful spectroscopic 
follow-up campaign and (b) the survey selection function 
and implied Lya nebula number density are presented in 
the two subsequent papers in this series (Prescott et al. 
2012b and Prescott et al. 2013, in preparation; hereafter. 
Paper II and III). 

We assume the standard ACDM cosmology {CIm = 0.3, 
ft A = 0.7, h = 0.7); V' corresponds to a physical scales 
of 8.5-7.6 kpc for redshifts of z ^ 2 — 3. All magnitudes 
are in the AB system. 

2. SURVEY DESIGN 

While most of the early examples of the Lya neb- 
ula class were found via narrow-band surveys of known 
galaxy over densities, one of the largest Lya nebulae 



2 



Prescott et al. 



was discovered in a very different manner (LABdOS; 
iDev et aI]l2QQ5f ). This Lya nebula at z ~ 2.7 came to 
light during a study of strong 24/im sources detected by 
the Spitzer Space Telescope in the NDWFS Bootes field. 
Strikingly, this source had a diffuse, extended morphol- 
ogy and very blue colors in deep broad-band imaging. 
The authors suspected and spectroscopically confirmed 
that Lya line and blue continuum emission were domi- 
nating the broad-band flux and that the diffuse morphol- 
ogy indicated the presence of a giant, spatially extended 
Lya nebula. 

Inspired by this discovery, we developed a search for 
Lya nebulae designed to work on deep broad-band imag- 
ing and applied it to the ~9 square degree NDWFS 
Bootes field dataset. Using a deep, wide-area, ground- 
based imaging survey like NDWFS allows us to probe an 
enormous comoving volume with existing data and signif- 
icantly reduce the amount of new telescope time required 
to complete the program. In that sense, our survey is 
complementary to the more sensitive but smaller volume 
narro w-band Lya neb ulae surveys (^e.g.. iMatsuda et al. 
[2004';"Saito et al."2006"; "Smith & Jaxvis|(2007|; Kang: et al. 
^009, 2010; Matsuda et al. 20Ul). In this section we de- 
scribe the NDWFS dataset and introduce the search algo- 
rithm used to select a sample of Lya nebulae candidates. 

2.1. NDWFS Broad-band Data 

Our survey technique was built and tested using the 
broad-band optical imaging from NDWFS, data which 
are available through the NO AO Science Archive Q The 
Bw and i?-band images of the 9.4 square degree Bootes 
field were obtained using the Mayall 4m Telescope and 
the MOSAIC 1 prime focus camera. The Bootes field 
was originally selected as a field suitable for a deep extra- 
galactic survey because of its low 100/im background and 
N(IIl) column density along the line of sight, and for be- 
ing well positioned from Kitt Peak National Observatory 
during a period of the year with historically good observ- 
ing conditions. The entire Bootes field was covered using 
stacked observations of 27 overlapping MOSAIC point- 
ings. The resulting median seeing was 1.10'' and 1.10'' 
in the Bw and R bands, respectively, and the median 
point-source depths were ^26.3 and 25.2 AB mag (5<j; 2" 
diameter apertures). The median la Bw surface bright- 
ness limit is 28.9 mag arcsec"^ (1.1" diameter apertures), 
which corresponds (in the case of pure line emission) to 
a line flux surface brightness limit of 2.48 x 10~^^ erg s~^ 
cm~^ arcsec"^. The data were reduced by the NDWFS 
team using 

2.2. The Search Algorithm 

The goal of our systematic survey was to find large, 
luminous Lya nebulae similar to LABd05. The search 
pipeline uses a blue broad-band image — in the case 
of the NDWFS, the 5v\/-band — as the primary search 
image and is tuned to find diffuse, spatially extended, 
low surface brightness objects. It is most sensitive to 
extended sources for which a bright Lya emission line 
boosts the broad-band flux relative to the very dark sky 
in the blue (Figure [1]), i.e., the largest and brightest Lya 



nebulae. A redder band — the i?-band in the case of 
NDWFS — is used to derive secondary color information. 

Given the depth and the size of the imaging survey, 
in designing our search algorithm it was important to 
minimize the number of contaminant objects while re- 
taining objects likely to be Lya nebulae. Our search 
strategy therefore involved first subtracting off bright 
galaxies and compact objects and then using a wavelet 
deconvolution algorithm to select diffuse, spatially ex- 
tended, low surface brightness objects with blue colors. 
The search pipeline steps are described in detail below 
and illustrated using a flowcharlQ in Figure [2j Figures [3] 
and m show the pipeline at work on the image region near 
LABd05. 

2.2.1. Background and Halo Subtraction 

Step 1: The first step is to remove the sky back- 
ground using a backgrou nd map generated by S ource 
Extractor (Version 2.4.4: iBertin fc ArnoutsI [19961 ) . We 
used a high detection threshold (detect.thresh =10(7 
per pixel), set background parameters (back.size =64 
and BACK_FILTERSIZE =3), and required a minimum of 
5 connected pixels (detect.minarea =5). This back- 
ground map, which includes a contribution from the ha- 
los around bright stars, was subtracted from the search 
image. 

2.2.2. Bright Object Removal 

In order to remove bright stars and galaxies that may 
contaminate the final candidate list, we employed the 
following procedure. 

Step 2: We first ran Source Extractor with a high de- 
tection threshold (DETECT_THRESH=10cr per pixel) and 
a minimum of 5 connected pixels (detect_minarea=5) 
in order to generate a catalog of bright sources. 

Step 3: We then ran Source Extractor with a very 
low threshold (detect_thresh= Icr per pixel) in or- 
der to generate a deep segmentation map, i.e., a map of 
all object pixels, that would include pixels in the faint 
wings of bright stars and galaxies. This iteration used 
the ASSOC mode in which the catalog from Step 2 was 
used as an input to Source Extractor and only those de- 
tections that match were included in the catalog. De- 
tected sources were matched against the input catalog 
in a weighted fashion based on measured pixel position 
(x_iMAGE, y_image), the object flux (flux_auto), and 
a search radius (assoc_radius=2 pixels). Using the 
deep segmentation maps from this step and the image 
rms maps, we replaced regions identified with bright stars 
and galaxies with patches of sky noise. 

Step 4' We ran Source Extractor on the "cleaned" 
image from Step 3 but this time with an intermediate 
threshold (DETECT_THRESH=6cr per pixel, determined 
based on visual inspection) to further remove stars and 
galaxies. 

Step 5: We ran Source Extractor on the image from 
Step 4 with the ASSOC feature in order to generate the 
corresponding deep segmentation map. Again, the de- 
tections are replaced with sky noise. 

2.2.3. Faint Compact Object Removal 



^ NOAO Science Archive: http://archive.noao.edu/nsa/ ^ Flowchart created using Omnigraffle: 

^ NDWFS Data Processing: http: / /www. noao.edu/n oao/noaodeep/] |http: / / www.o mnigroup.com / products / omnigraffle / 
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Step 6: We next used unsharp masking to remove 
fainter objects with compact morphology. An unsharp 
masked image was produced by smoothing the image 
from Step 5 using an 11 x 11 pixel kernel and then sub- 
tracting the smoothed image from the original Step 5 
image. In this process, diffuse emission was removed, 
leaving only compact sources in the unsharp masked im- 
age. 

Step 7: We performed another ASSOC Source Extrac- 
tor run on the unsharp masked image in which we only 
allowed matched detections to be included in the catalog, 
i.e., anything detected in both the Step 5 image and the 
unsharp masked image was reported. We removed these 
object regions by replacing them with sky noise, and in 
this way removed additional faint compact sources from 
the image. 

2.2.4. Final Image Cleaning 

Step 8: To avoid artifacts (e.g., saturated columns, 
poorly-subtracted halos) from bright stars being selected 
as sources later in the pipeline, we replaced the regions 
around bright stars {B < 16 mag in the USNO-A2.0 
catalog) out to a radius of 150 pixels (~ 39^') with sky 
noise. This step removed ?^0.68 square degrees from the 
final survey area. 

2.2.5. Identification of Spatially Extended Objects 

Step 9: We then applied a wavelet decomposition al- 
gorithm to decompose the image into maps showing the 
power on a number of different spatial scales. 

The basic idea behind wavelet decomposition is to de- 
compose the image into a set of "wavelet planes," each 
of which contains the flux of sources with power on that 
scale and the sum of which is the original image (ignor- 
ing losses). The power of the wavelet decomposition ap- 
proach is to reduce the interference of small scale objects 
when looking for large scale structures. 

We use the wavelet decomposition code wvdecom^ to 
filter the image across six size scales {scalemin = 1, 
scalemax = 6, threshold = 3, thresholdmin = 2, 
iter = 10). Here we give a brief description of the al- 
gorithm; a detailed discussion of the code is given in the 
ZHTOOLS documentation. The wvdecomp algorithm uses 
an d trous wavelet kernel which is defined as the dif- 
ference between two functions (/^ and /i+i) that each 
roughly resemble a Gaussian of width 2*~^. On each 
scale, the kernel is sensitive to objects of the same size 
(« 2-1). 

The algorithm begins by convolving the image with the 
smallest wavelet kernel. Any features that are insignifi- 
cant in the convolved image are zeroed and the result is 
subtracted from the original image, causing a majority 
of the flux from small features comparable in size to the 
smallest wavelet kernel to be removed. The process is 
repeated with larger and larger wavelet kernels, and the 
results from each step are saved as wavelet power maps 
or "wavelet planes." To minimize the loss of information 
at each step, two thresholds are specified in running the 
code: a detection threshold, imax^ and a filtering thresh- 
old, 

^mini both defined in terms of the image noise. A 

^ wvdecomp was written by Alexey Vikhlinin; 
ZHTOOLS documentation and source code: 

http: / / hea-www.harvard.edu/RD / zhtools / 



source detected on a given scale above irnin is counted as 
significant only if the source local maximum also exceeds 

Step 10: We used the known Lya nebula LABd05 to 
empirically determine the appropriate scale for large Lya 
nebulae. Using one final run of Sour ceExtr actor, we gen- 
erate a catalog of sources in the wavelet power map and 
select those with wavelet power peaks above 4cr, chosen 
to minimize the number of total candidates while pre- 
serving sources like LABd05. 

2.2.6. Final Catalog Cleaning and Candidate Prioritization 

Step 11: Sources drawn from low exposure time regions 
of the survey (<4000 sec) as well as those found in regions 
that had been flagged in previous steps were removed 
from the candidate list. This step removed ^0.28 square 
degrees from the final survey area. 

Step 12: Postage stamps of candidates were inspected 
visually and classified by eye using morphological flags, 
where STAR denotes a star halo, galaxy denotes a 
galaxy halo, tidal/arm indicates tidal tails or spiral 
arms, GROUP corresponds to a tight grouping of compact 
sources selected as a single large source, and diffuse in- 
dicates spatially extended diffuse emission. The 2266 
sources selected by the morphological search pipeline 
were categorized as follows: 738 STAR sources, 491 
GALAXY sources, 170 tidal/arm sources, 459 group 
sources, and 408 diffuse sources. We removed those 
objects categorized as STAR, galaxy, and tidal/ ARM 
from the candidate sample (1399 sources), since these 
are easy to classify unambiguously by eye. As the dis- 
tinction between categories GROUP and diffuse is more 
subjective, both category GROUP and diffuse objects 
were included in the final sample (867 sources). 

2.2.7. Candidate Selection 

Step 13: The final sample was then prioritized based 
on Bw — R color and size, giving priority to blue can- 
didates with large sizes. Figure [5] shows the size and 
color distribution of the candidate sample. We measured 
Bw — R colors in the original images using large aper- 
tures (30 pixel, 7.7'' diameter, i.e., roughly 60 kpc; aper- 
ture sky subtraction). The size of each source, for the 
purposes of selecting candidates, was taken to be the size 
measured in the wavelet power map (isoarea.world in 
square degrees), as determined by SourceExtractor dur- 
ing the final search {Step 10). While representative, this 
"wavelet size" should not be taken to be the true size 
of the object. The wavelet size typically underestimates 
the nebular extent in cases where some portion of that 
object has been rejected during the pipeline (due to com- 
pact sources within or adjacent to the diffuse emission). 
In addition, as we are attempting to detect line emission 
within the broad Bw band, the measured Bw sizes will 
underestimate the true nebular size. We investigate this 
in more detail in Paper II. 

The majority of the morphologically-selected candi- 
dates are red {Bw — R ^ 0.5), but a small subset show 
blue colors, which could either result from a strong emis- 
sion line within the Bw band (Lya at 1.9 < z < 2.9), 
from very blue continuum emission, or a combination. 
Figure [6] shows the expected Bw — R color for a model 
source with a given redshift and Lya equivalent width 
generated by adding a Gaussian emission line to a flat 
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spectrum source {fx ex A~^). Unlike narrow-band Lya 
nebula surveys, our survey is not formally equivalent 
width limited. Other than at the high redshift end of the 
survey {z > 2.7), sources of any equivalent width will fall 
within our color selection window as long as they are spa- 
tially extended enough and bright enough in Bw to be 
selected via the morphological search pipeline. However, 
given the distribution of field galaxy colors we expect 
that sources with bluer colors are more likely to be Lya 
nebulae. Given the number of candidates we expected 
to be able to target during our spectroscopic follow-up 
campaign, we therefore defined our first priority region 
as Bw - R < 0.45 and isoarea_world> 34 arcsec^ 
and a second priority region as 0.45 < Bw — R < 0.65 
and Log(isOAREA_wORLD)> 0.934 x {Bw - i^) + 1.11 
arcsec^ (Figure [5j). Our spectroscopic follow-up observa- 
tions (using the MMT/Blue Channel Spectrograph with 
a full wavelength range of AA = 3100-8320A) are able to 
detect line emission from Lya nebulae down to z ^ 1.6, 
leading to a total survey redshift range ofl.6<z<2.9. 
The details of our selection function will be treated in 
Paper III. 

3. CANDIDATE SAMPLE 

Our final candidate sample consists of 39 first priority 
and 40 second priority sources over a search area of 8.5 
square degrees (Table [T]). Of the 39 first priority sources, 
32 have a morphological category of GROUP and 7 are in 
the category DIFFUSE; for the 40 second priority sources, 
31 are in the GROUP category and 9 in the DIFFUSE cate- 
gory. The distribution of Lya nebulae candidates on the 
sky is shown in Figure [71 and postage stamps of the Bw^ 
i?, and /-band imaging of all candidates are shown in 
Appendix [Al Due to the fact that we included both cate- 
gory GROUP and DIFFUSE sources in the final catalog, the 
candidates span a range of morphologies from clumpy to 
diffuse. In particular, a number of the very blue sources 
appear to be a close grouping of compact sources selected 
as one. While at first glance these compact groups do not 
have the diffuse appearance we expect for Lya nebulae, it 
could be unwise to remove them because the full morpho- 
logical range spanned by the Lya nebula class is not well- 
known and some Lya nebulae are known to include com- 
pact sources (e.g., young galaxies; Prescott et al. 2012c, 
submitted). For these reasons we did not apply any fur- 
ther morphological selection criteria to the main sample. 
In addition to our 79 systematically-selected first and 
second priority samples, we also flagged 6 additional can- 
didates with promising morphologies from outside these 
regions as a test of our color-selection approach. Of the 
third priority sources, which are included for reference in 
Table [H and Appendix [Al 5 have a morphological cate- 
gory of DIFFUSE and 1 is in the category GROUP. 

Based on previous narrow-band Lya nebulae surveys, 
we can predict the number of Lya nebulae we would 
expect to find within our full survey volume. Figure [8] 
shows the cumulative number of Lya nebulae expected 
assuming a constant comoving volume density as a func- 
tion of redshift of 0.5 x 10"^ - 3. 6 x 10"^ M pc"^ (the 
range spanned by the estimates of lYang et al.l l2009) and 
a 100% detection success rate. This non-evolving model 
would predict approximately 60-400 Lya nebulae over 
the redshift range of our survey (1.6 < z < 2.9). As 
the comoving volume density of Lya nebulae is unlikely 



to remain constant over cosmic time (i Keel et al.l [20091 ). 
for comparison we also plot the cumulative number of 
Lya nebulae expected within our survey area assuming 
that the comoving volume density increases linearly as a 
function of redshift (dashed lines). In reality the selec- 
tion function, completeness, and efliciency of our survey 
are differen t from , that of narrow-band surveys such as 
[Yang et"al] (|2009l ) and the evolution of the Lya nebula 
comoving volume density is unknown; we address these 
issues in more detail in Paper III. 

3.1. Potential Contaminants 

Due to the wide bandpass of the Bw filter, there 
are several populations that may contaminate our Lya 
nebula candidate sample. The morphological selection 
method will tend to pick out any source of low surface 
brightness extended emission that cannot be easily iden- 
tified, e.g., a low surface brightness galaxy (LSB) or low 
surface brightness nebulae within the Galaxy. In princi- 
ple, lower redshift sources with [Oll]AA3727,3729 emis- 
sion in the selection band {Bw) are a contaminant popu- 
lation, much as they are in Lya-emitting galaxy surveys. 
In addition, galaxies and Lya nebulae in the redshift 
desert (i.e., 1.2 < z < 1.6) may be selected as candi- 
dates in our survey if they exhibit extended blue contin- 
uum emission or a close grouping of compact sources. In 
these cases, Lya would be positioned blueward and [Oil] 
redward of the wavelength range of our follow-up spec- 
troscopy, leading us to detect faint, spatially extended, 
blue continuum emission with no emission lines to pro- 
vide confirmation of the redshift. 

While this range of different populations could contam- 
inate our sample in principle, not all of them will con- 
tribute substantially in practice. We do not expect to 
find many LSBs because the Bw — R colors of Lya neb- 
ulae are likely to be bluer than typical LSB colors (i.e., 
Bw - R=IA9± 0.60: Haberzettl et al. 2007): any blue 
LSBs would be expected to show optical emission lines 
and therefore easily identified in follow-up spectroscopy. 
If compact groupings of [Oil] -emitters or extended [Oil] 
nebulae were to be selected via our morphological search, 
the position of [Oil] in the Bw, [Om] the R band, and 
Ha in the R or I band would give the object redder col- 
ors, which would not be included in the final sample. 
On the other hand, a primary contaminant population 
will be galaxies or Lya nebulae in the redshift range 
1.2 < z < 1.6 that appear in optical spectroscopy as spa- 
tially extended continuum-only sources. From Figure [HI 
we would expect of order 20-180 such sources, or roughly 
25% of the predicted Lya nebulae sample over the en- 
tire redshift range (1.2 < 2; < 2.9), assuming no evolu- 
tion in the intrinsic number density. With line emission 
outside the optical window, however, follow-up observa- 
tions with HST will be necessary to confirm whether or 
not any continuum-only sources are in fact lower redshift 
Lya nebulae. 

4. DISCUSSION 

4.1. A Successful Broad-band Lya Nebula Survey 

Our systematic broad-band Lya nebula search was de- 
signed to find luminous Lya nebulae similar to LABd05. 
Subsequent spectroscopic follow-up of a subset of our 
morphological candidate sample showed that, in addi- 
tion to recovering the previously-discovered large Lya 
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nebula z ^ 2.66 fLABdQ5:l Dev et alJ[2QQ5l ). this novel 
approach was able to successfully locate 4 new large 
hya nebulae over the NDWFS Bootes field z ^ 1.67, 
z ^ 2.27, z ^ 2.14, and z ^ 1.88 (PRGl, PRG2, 
PRG3, PRG4; |Prescptt"2009) . Of particular note, our 
broad-band search identifie d the first spat ially extended 
Lya+Hell nebula fPRGl: iPrescott et al.i r2QQ9): located 
at z ^ 1.67, this system is the lowest redshift Lya neb- 
ula known, and the strong Hell emission suggests that it 
may contain low metallicity gas. The 5 confirmed Lya 
nebulae are shown for reference in Figures [5]l6l [71 and 
[9l A full discussion of the spectroscopic follow-up cam- 
paign is presented in Paper II, and the selection function 
and implied space density for this survey are discussed 
in Paper III. 

4.2. Comparisons to Other Surveys 

While broad-band imaging has been used in recent 
years to isolate samples of Ly g- emittin g and Lya- 
absorbing galaxies based on colors (jCooket 2009) and to 
find strong line-emitting nebulae at low redsh i fts through 
the G alaxy Zoo project (jLintott et al.l 120091 : iKeel et aTl 
our search is the first time broad-band data has 
been used as the basis for a semi-automated system- 
atic morphological survey for high redshift Lya nebu- 
lae. Our Lya nebula survey is also fundamentally differ- 
ent and complementary to standard narrow-band surveys 
for Lya nebulae (Table [2]). The pilot survey of the ND- 
WFS Bootes field has shown that it is possible to find 
large Lya nebulae using deep broad-band imaging data. 
There are many advantages to this unusual approach. 
A broad-band search is able to efficiently cover enor- 
mous comoving volumes that would be observationally 
expensive when using a narrow-band filter. A narrow- 
band survey using a 100 A filter at z 2.3 would need to 
search over 100 square degrees to reach the same comov- 
ing volume covered in our broad-band survey of the ND- 
WFS Bootes field 10^ hfg Mpc"^). In addition, our 
technique builds on existing datasets, which means that 
our search algorithm can be easily modified to find Lya 
nebulae within any deep, broad-band, wide-field imaging 
survey, leading to substantially lower observational over- 
head. Finally, while the early, successful Lya surveys 
often targeted known overdensities and were limited to 
smaller areas (iSteidel et aDjgOOOt iMatsuda et al.l 120041 : 
iFrancis et al.l 1996; Palun as et al l [200^ . introducing a 
bias to the number density estimates, our 8.5 square de- 
gree survey, by contrast, is unbiased in terms of envi- 
ronment, spanning a Az ~ 1.3, a comoving transverse 
distance of ^290 h^Q^ Mpc (at the redshift midpoint, 
z ^ 2.25), and a comoving line-of-sight distance of ^1700 
hjQ Mpc. 

At the same time, a drawback of using broad-band 
data to select line-emitting sources is the necessity of re- 
stricting our search to the largest, most luminous Lya 
nebulae and of removing compact sources in favor of dif- 
fuse emission, in order to reduce the number of contam- 
inants. Thus, while our search successfully recovered a 
sample of Lya nebulae, a handful of other (fainter or 
more compact) Lya nebulae that were identified by tra- 
ditional narrow-band searches in Bootes were not recov- 
ered by our broad-band search (Table [3]). 

In Figure [9l we show the Bw isophotal area and mean 



surface brightness estimates for the candidate sample, as 
derived using SourceExtractor (detect_thresh= 28.9 
mag arcsec"^, detect_minarea= 5) on the original Bw 
imaging, for the Lya nebulae confirmed by our survey 
(Paper II), and for other known Lya nebula sources in 
the field. The sources our survey did not recover fall 
into two categories: either the Lya emission surrounds 
a much brighter compact continuum source ( "core-halo" 
morphology) or the diffuse Bw emission is simply too 
faint to be selected by our search algorithm. In the first 
category are the four Lya nebulae found bv I Yang et al.l 
(2009) using a traditional narrow-band survey of ND- 
WFS Bootes; all four consist of bright continuum sources 
(AGN in two cases) surrounded by Lya halos that are 
undetected in the broad-band Bw imaging. The bright 
compact central sources caused these objects to be re- 
jected and replaced with sky noise during Steps 2-3 or 
Steps 6-7 of the search pipeline. Intermediat e-band imag- 
ing with the Subaru (jPrescott et al.l I2008D and Mayah 
Telescopes revealed three additional Lya nebulae that 
were not detected using the broad-band search algo- 
rithm. One case (PI) was too faint to be selected by our 
search algorithm, and two other cases (P2 & P3) showed 
core-halo morphologies with confirmed Lya emission (at 
z ^ 2.6 and z ~ 1.9) in the midst of several bright contin- 
uum sources. With the compact sources removed by the 
search pipeline, the wings of the Lya emission were still 
detected in the Bw imaging in each case, but at weaker 
wavelet power limits than probed by the current survey. 

As our search was designed to select large, luminous, 
and diffuse Lya nebulae, it is not surprising that more 
compact or fainter Lya nebulae were not selected. In fu- 
ture work, it will be important to understand this range 
of Lya nebula morphologies and whether they reflect 
fundamentally different sources or simply different evo- 
lutionary phases of the same underlying phenomenon. 

5. CONCLUSIONS 

We have designed an innovative systematic search for 
large Lya nebulae using deep broad-band data. Our 
technique is designed to find the largest and brightest 
Lya nebulae and is able to probe enormous comoving 
volumes 10^ h^Q Mpc^) using existing deep broad- 
band datasets. In this paper, we presented details on the 
survey design and our sample of Lya nebula candidates 
selected from the NDWFS Bootes Field. Paper II of this 
series presents the spectroscopic follow-up and the con- 
firmed Lya nebula sample, and Paper III will address 
the details of the survey selection function and the im- 
plications for the space density of Lya nebulae. In the 
future, we will build on the search algorithm developed 
in this work in order to find Lya nebulae in other deep 
and wide imaging fields and increase the sample of known 
Lya nebulae available for detailed follow-up study. 

The authors are grateful to Alexey Vikhlinin for the 
use of the code wvdecomp. This research draws upon data 
from the NDWFS as distributed by the NOAO Science 
Archive. NOAO is operated by the Association of Uni- 
versities for Research in Astronomy (AURA), Inc. under 
a cooperative agreement with the National Science Foun- 
dation. M. P. was supported by a NSF Graduate Re- 
search Fellowship and a TABASGO Prize Postdoctoral 
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Fellowship. A. D. and B. T. J.'s research is supported by 
NO AO, which is operated by AURA under a cooperative 
agreement with the National Science Foundation. 
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Table 1 

Lya Nebula Candidates 





Candidate Name 


Right Ascension 


Declination 


Wavelet Size 


Morphological 




Bw -R^ 


Priority^ 








(J2000 hrs) 


(J2000 deg) 


(arcsec^) 


Flag^ 


(AB mag) 


(AB mag) 




1) 


NDWFS 


J143004.7+353509 


14:30:04.687 


35:35:09.06 


134.0 


DIFFUSE 


23.79+0.05 


0.61 


2d 


2) 


NDWFS 


J143006.9+353437 


14:30:06.864 


35:34:36.73 


108.2 


DIFFUSE 


23.05+0.03 


0.80 


3d 


3) 


NDWFS 


J142846.2+330819 


14:28:46.228 


33:08:19.42 


89.9 


DIFFUSE 


23.00+0.03 


0.89 


3 


4) 


NDWFS 


J143436.4+330406 


14:34:36.448 


33:04:05.62 


72.1 


GROUP 


22.90+0.02 


0.53 


2 


5) 


NDWFS 


J142855.0+353022 


14:28:54.998 


35:30:21.70 


64.9 


DIFFUSE 


23.26+0.03 


0.52 


2 


6) 


NDWFS 


J143153.1+351436 


14:31:53.066 


35:14:36.45 


62.6 


GROUP 


23.31+0.04 


0.53 


2 


7) 


NDWFS 


J143055.2+352855 


14:30:55.192 


35:28:54.76 


61.6 


GROUP 


23.30+0.04 


0.43 


1 


8) 


NDWFS 


J143057.7+354547 


14:30:57.662 


35:45:47.26 


59.6 


GROUP 


23.09+0.03 


0.43 


1 


9) 


NDWFS 


J143842.0+335325 


14:38:41.956 


33:53:24.79 


59.5 


GROUP 


23.14+0.04 


0.63 


2 


10) 


NDWFS 


J143410.9+331731 


14:34:10.948 


33:17:30.80 


59.1 


DIFFUSE 


23.18+0.03 


0.42 


1^ 


11) 


NDWFS 


J142729. 9+341125 


14:27:29.894 


34:11:25.22 


56.2 


GROUP 


23.48+0.05 


0.48 


2 


12) 


NDWFS 


J143551. 3+333421 


14:35:51.278 


33:34:21.03 


55.9 


GROUP 


23.38+0.05 


0.51 


2 


13) 


NDWFS 


J142909.3+342542 


14:29:09.278 


34:25:42.38 


55.8 


GROUP 


23.41+0.04 


0.58 


2 


14) 


NDWFS 


J143512.3+351109 


14:35:12.336 


35:11:08.62 


55.7 


DIFFUSE 


23.53+0.04 


0.56 


2f 


15) 


NDWFS 


J142906.3+341130 


14:29:06.271 


34:11:30.26 


53.5 


GROUP 


23.47+0.04 


0.48 


2 


16) 


NDWFS 


J142856.0+331937 


14:28:56.013 


33:19:36.51 


53.5 


GROUP 


23.38+0.04 


0.60 


2 


17) 


NDWFS 


J143511. 0+335543 


14:35:11.030 


33:55:43.10 


53.3 


GROUP 


23.12+0.03 


0.53 


2 


18) 


NDWFS 


J143222. 8+324943 


14:32:22.768 


32:49:42.67 


53.1 


DIFFUSE 


23.35+0.04 


0.53 


2 


19) 


NDWFS 


J143302.8+354729 


14:33:02.844 


35:47:29.43 


52.7 


GROUP 


23.42+0.04 


0.50 


2 


20) 


NDWFS 


J143641.4+341510 


14:36:41.373 


34:15:10.26 


52.1 


GROUP 


23.50+0.04 


0.55 


2 


21) 


NDWFS 


J143230.2+334105 


14:32:30.244 


33:41:05.46 


50.5 


GROUP 


23.25+0.03 


0.56 


2 


22) 


NDWFS 


J142819.8+344657 


14:28:19.840 


34:46:57.14 


50.4 


GROUP 


23.16+0.03 


0.61 


2 


23) 


NDWFS 


J143115. 1+354148 


14:31:15.055 


35:41:47.76 


50.1 


GROUP 


23.54+0.05 


0.55 


2 


24) 


NDWFS 


J142614.7+344434 


14:26:14.714 


34:44:34.22 


49.9 


GROUP 


23.51+0.05 


0.43 


1 


25) 


NDWFS 


J142929.2+354711 


14:29:29.198 


35:47:10.78 


49.5 


GROUP 


23.17+0.03 


0.52 


2 


26) 


NDWFS 


J142622.9+351422 


14:26:22.905 


35:14:22.02 


49.5 


DIFFUSE 


23.20+0.04 


-0.49 


IS 


27) 


NDWFS 


J143055.7+340502 


14:30:55.713 


34:05:01.71 


49.3 


DIFFUSE 


23.04+0.03 


0.57 


2 


28) 


NDWFS 


J143236.9+351800 


14:32:36.883 


35:18:00.32 


49.1 


GROUP 


23.42+0.06 


0.57 


2 


29) 


NDWFS 


J142526.3+335112 


14:25:26.332 


33:51:12.16 


48.7 


GROUP 


23.35+0.04 


0.40 


1 


30) 


NDWFS 


J143407.5+334141 


14:34:07.476 


33:41:40.56 


48.3 


GROUP 


23.43+0.05 


0.26 


1 


31) 


NDWFS 


J142547. 1+334454 


14:25:47.126 


33:44:54.13 


48.1 


DIFFUSE 


22.98+0.03 


0.45 


1 


32) 


NDWFS 


J143459.7+333749 


14:34:59.702 


33:37:48.93 


48.1 


GROUP 


23.68+0.06 


0.48 


2 


33) 


NDWFS 


J142714.8+343155 


14:27:14.791 


34:31:54.55 


47.7 


GROUP 


23.38+0.04 


0.52 


2 


34) 


NDWFS 


J143128.2+352658 


14:31:28.245 


35:26:57.91 


47.3 


DIFFUSE 


23.06+0.03 


0.56 


2 


35) 


NDWFS 


J142906.4+352432 


14:29:06.398 


35:24:32.22 


47.2 


GROUP 


23.50+0.05 


0.54 


2 
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Table 1 — Continued 





Candidate Name 


Right Ascension 


Declination 


Wavelet Size 


Morphological 






Bw - 


Priority'^ 








(J2000 hrs) 


(J2000 deg) 


(arcsec^) 


Flag^ 


(AB ma 


%) 


(AB mag) 




36) 


IN U VV r o 


J i4o4zz.o-t-o0i0o4 


14:34:22.348 


35:15:34.27 


46.3 


GROUP 


23.50+0, 


.04 


0.51 


2 


37) 


IN U VV r o 


Ti/iqq-|n QiQQC^nP^T 


14:33:10.346 


33:50:56.94 


45.9 


GROUP 


23.82+0, 


.07 


0.58 


2 


38) 


IN U VV r o 


T1/1Q[^OQ Q_I_QC^/1QQ/1 

J i4oozo.y-t-oo4yo4 


14:35:23.856 


35:49:34.39 


45.0 


GROUP 


23.93+0, 


.07 


0.42 


1 


39) 


IN U VV r o 


J i4oU0o.0H-o0zUU / 


14:30:53.488 


35:20:06.68 


45.0 


GROUP 


23.02+0, 


.04 


0.43 


1 


40) 


NDWl S 


J 142653.2+343855 


14:26:53.172 


34:38:55.39 


44.8 


GROUP 


23.53+0, 


.06 


-0.67 


1^ 


41) 


NDWFS 


J142758. 6+354428 


14:27:58.617 


35:44:28.42 


44.4 


GROUP 


23.42+0, 


.03 


0.50 


2 


42) 


JNDWl S 


J 142634. 1+334954 


14:26:34.106 


33:49:53.61 


44.3 


GROUP 


23.54+0, 


.04 


0.51 


2 


43) 


JNDWl S 


J 143435.9+350710 


14:34:35.932 


35:07:09.91 


44.2 


GROUP 


23.72+0, 


.05 


0.57 


2 


44) 


AT TAT H TTTIO 


J 142927. 8+34590O 


14:29:27.837 


34:59:06.14 


43.1 


GROUP 


23.18+0, 


.03 


1.18 


3 


45) 


NDWl S 


J 142745.2+352949 


14:27:45.194 


35:29:48.98 


42.9 


GROUP 


23.48+0, 


.04 


0.23 


1 


46) 


NDWl S 


J 142826.0+352116 


14:28:25.999 


35:21:15.62 


42.7 


GROUP 


23.58+0, 


.04 


0.48 


2 


47) 


NDWFS 


J142535. 2+324934 


14:25:35.205 


32:49:34.17 


42.6 


GROUP 


23.54+0, 


.05 


0.49 


2 


48) 


N U W 1^ D 


J 142900.0+050927 


14:29:59.978 


35:39:26.89 


42.4 


GROUP 


23.36+0, 


.03 


0.47 


2 


49) 


IN U VV r O 


J i4ZoZi.o+o0zUUi 


14:28:21.259 


35:20:00.92 


42.4 


GROUP 


23.42+0, 


.04 


0.25 


1 


50) 


IN U VV r o 


J ±404u±.u-|-ooizzy 


14:34:01.034 


35:12:29.12 


42.0 


GROUP 


23.91+0, 


.06 


0.19 


1 


51) 


IN U VV r o 


J i4oo4o.y-t-o04i0 I 


14:33:48.885 


35:41:56.68 


41.7 


GROUP 


23.74+0, 


.05 


-0.08 


1 


52) 


IN U VV r o 


J i4o < UO.O-t-OOOOOo 


14:37:06.588 


33:56:52.65 


41.3 


DIFFUSE 


23.22+0, 


.03 


0.52 


2 


53) 


IN U VV r o 


J i4z440.U+o04 / 14 


14:24:45.998 


35:47:14.06 


40.5 


GROUP 


23.68+0, 


.05 


0.19 


1 


54) 


IN U VV r o 


J ±4zOzZ.O-|-OZ04z4 


14:25:22.339 


32:54:23.65 


39.9 


GROUP 


23.40+0, 


.04 


0.48 


2 


55) 


IN U VV r o 


J ±40o00.Z-|-04U00y 


14:38:33.249 


34:03:59.32 


39.6 


GROUP 


23.74+0, 


.06 


0.52 


2 


56) 


IN U VV r o 


J i40±4±.0-|-O0z±iU 


14:31:41.253 


35:21:09.86 


39.5 


DIFFUSE 


23.17+0, 


.05 


0.48 


2 


57) 


IN U VV r o 


J i404±0.4-|-000U0o 


14:34:16.420 


33:30:58.10 


39.3 


GROUP 


23.26+0, 


.05 


0.46 


2 


58) 


IN U VV r o 


J i4z0iD.D-|-0Z4000 


14:25:16.629 


32:43:35.47 


39.0 


GROUP 


23.13+0, 


.04 


0.14 


1 


59) 


NDWl S 


J 143412.7+332939 


14:34:12.722 


33:29:39.19 


38.7 


DIFFUSE 


23.40+0, 


.07 


0.22 


ih 


60) 


NDWl S 


J 142441.9+332532 


14:24:41.940 


33:25:31.94 


38.4 


GROUP 


23.52+0, 


.04 


0.38 


1 


61) 


NDWl S 


J 143438.7+350420 


14:34:38.700 


35:04:19.84 


38.4 


DIFFUSE 


23.72+0, 


.05 


0.46 


2 


62) 


AT TAT H TTnO 


J 142832. 7+O54105 


14:28:32.733 


35:41:04.56 


37.8 


GROUP 


23.72+0, 


.05 


-0.23 


1 


63) 


AT TAT H TTTIO 


J 14ooo5.8+o4450o 


14:33:35.764 


34:45:02.84 


37.5 


GROUP 


23.92+0, 


.06 


0.33 


1 


64) 


ATTAT H 7TnO 

NDWl^ b 


J 143203.8+351855 


14:32:03.760 


35:18:54.86 


37.4 


GROUP 


23.46+0, 


.05 


0.46 


2 


65) 


NDWFS 


J143207. 2+343101 


14:32:07.224 


34:31:01.34 


37.3 


DIFFUSE 


23.30+0, 


.03 


0.81 


3 


66) 


NDWFS 


J142539. 9+344959 


14:25:39.859 


34:49:59.19 


37.1 


GROUP 


23.57+0, 


.04 


0.14 


1 


67) 


NDWl S 


J 142819.7+325449 


14:28:19.735 


32:54:48.70 


37.1 


GROUP 


23.44+0, 


.03 


0.44 


1 


68) 


NDWl S 


T"1 A C\ A A ^ f\ 1 O/I/IOO/I 

J 142443.9+344834 


14:24:43.869 


34:48:34.45 


36.9 


GROUP 


23.64+0, 


.05 


0.37 


1 


69) 


AT TAT T TTTIO 

NDWl b 


J 143202. D+O51904 


14:32:02.568 


35:19:04.22 


36.6 


DIFFUSE 


22.99+0, 


.04 


0.42 


1 


70) 


ATTAT T 7Tr^O 

NDWl b 


J 142753.8+341204 


14:27:53.762 


34:12:04.10 


36.5 


DIFFUSE 


23.29+0, 


.04 


0.45 


1 


71) 


TVTTATS TTPC 

N D W 1 b 


J 142600. 8+O50252 


14:26:00.842 


35:02:52.36 


36.1 


DIFFUSE 


23.87+0, 


.06 


0.63 


3 


72) 


NDWFS 


J142643. 9+340937 


14:26:43.850 


34:09:36.82 


36.1 


GROUP 


23.40+0, 


.03 


0.03 


1 


73) 


NDWl S 


J 142722.4+345225 


14:27:22.408 


34:52:24.74 


35.9 


GROUP 


23.71+0, 


.05 


0.29 


1 


74) 


NDWl S 


J 142620.0+340427 


14:26:19.982 


34:04:27.01 


35.7 


GROUP 


23.51+0, 


.04 


-0.75 


1 


75) 


NDWFS 




14:27:07.840 


34:47:48.84 


35.4 


DIFFUSE 


23.51+0, 


.05 


0.43 


1 


76) 


NDWFS 


J142710.4+324842 


14:27:10.447 


32:48:41.76 


35.4 


GROUP 


23.09+0, 


.04 


0.41 


1 


77) 


NDWFS 


J143232. 5+351534 


14:32:32.484 


35:15:34.23 


35.3 


GROUP 


23.45+0, 


.07 


0.30 


1 


78) 


NDWFS 


J143651. 3+342107 


14:36:51.312 


34:21:07.38 


35.1 


GROUP 


23.71+0, 


.05 


0.28 


1 


79) 


NDWFS 


J142746. 3+344544 


14:27:46.348 


34:45:44.02 


35.0 


GROUP 


23.39+0, 


.04 


0.34 


1 


80) 


NDWFS 


J142548. 3+322957 


14:25:48.283 


32:29:56.58 


35.0 


GROUP 


23.83+0, 


.06 


0.12 


1 


81) 


NDWFS 


J142735. 5+342332 


14:27:35.479 


34:23:32.38 


35.0 


GROUP 


23.55+0, 


.04 


-0.26 


1 


82) 


NDWFS 


J142449. 8+324743 


14:24:49.761 


32:47:42.61 


34.9 


GROUP 


24.09+0, 


.09 


-0.74 


1 


83) 


NDWFS 


J142732. 5+341213 


14:27:32.520 


34:12:13.39 


34.4 


GROUP 


23.62+0, 


.05 


-0.62 


1 


84) 


NDWFS 


J142802.8+350933 


14:28:02.803 


35:09:33.30 


34.2 


GROUP 


23.80+0, 


.06 


0.35 


1 


85) 


NDWFS 


J142533.0+343912 


14:25:32.966 


34:39:11.95 


28.8 


DIFFUSE 


23.67+0, 


.06 


0.26 


3 



^ Morphological flag discussed in Section 12.2.61 GROUP - tight 
grouping of compact sources, diffuse - spatially-extended diffuse 
emission. 

^ Magnitudes and colors were measured using large 30 pix {7.7") 
diameter apertures. 

^ Priority assignments are discussed in Section 12.2.71 
^ The two largest candidates are in fact part of a single diffuse, 
asymmetric, and mostly linear structure. The unusual morphology 
and large size suggest that this is more likely some sort of Galactic 
nebula rather than an extragalactic source. The optical spectrum 
(Paper II) does not show any obvious emission features, suggesting 
that this might be a reflection nebula. 



^ Our search recovered LABd05, the Lyo; nebula found previ- 
ously bv IDev et al.l (|2005l) . The listed coordinates are those re- 
turned by the search algorithm. 

^ Our search discovered PRGl, a Lyo; nebula at z ~ 1.67 
(Prescott et al. 2009, Paper II). 

s Our searc h discovered PRG2, a Lyo; nebula at z ~ 2.27 
(IPrescottil2009l . Paper II). 

^ Our search discovered PRG3, a Lyo; nebula at z ~ 2.14 
(jPrescott 2009, Paper II). 

^ Our search discovered PRG4, a Lyo; nebula at 2; ~ 1.88 
(IPrescottil2009l . Paper II). 
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Table 2 

Wide- Area Lya Nebula Surveys 



Survey 


Field^ 


Redshift 


Az 


Area 
(deg2) 


Comoving Volume 
(10^ hf(f Mpc^) 


Matsuda et al. 2004 


SSA22 


3.09 


0.06 


0.20 


0.14 


Palunas et al. 2004 


J2143-4423 


2.38 


0.04 


0.56 


0.30 


Saito et al 2006 


SXDS 


3.34 


0.20 


0.23 


0.53 




SXDS 


3.72 


0.22 


0.23 


0.56 




SXDS 


3.93 


0.24 


0.23 


0.60 




SXDS 


4.12 


0.24 


0.23 


0.59 




SXDS 


4.35 


0.26 


0.23 


0.62 




SXDS 


4.58 


0.28 


0.23 


0.65 




SXDS 


4.82 


0.26 


0.23 


0.59 


orQiin 66 jarvis zuu/ 


YVlVllVl-ijoo, ijii, or Ijo 


Z.oO 


U.Uo 


~±o 


/I 7 




AT" T\/n\/r TOd TTT CI T O 

XMM-LSS, LH, SFLS 


3.00 


0.14 


^.1 crb 


f^24.8 




XMM-LSS, LH, SILS 


3.12 


0.08 


?^15^ 


?^14.5 


Matsuda et al. 2009 


B3J2330+3927 


3.09 


0.06 


0.25 


0.18 


Ouchi et al. 2009 


SXDS 


6.56 


0.10 


1.00 


0.82 


Yang et al. 2009 


NDWFS Bootes 


2.30 


0.04 


4.82 


2.14 


Yane et al. 2010 


CDF-N 


2.30 


0.04 


0.28 


0.12 


CDF-S 


2.30 


0.04 


0.24 


0.11 




COSMOSl 


2.30 


0.04 


0.36 


0.16 




COSMOS2 


2.30 


0.04 


0.36 


0.16 


Matsuda et al. 2011 


SSA22 


3.09 


0.06 


1.06 


0.74 




SXDS 


3.09 


0.06 


0.60 


0.42 




GOODS-N 


3.09 


0.06 


0.24 


0.17 




SDF 


3.09 


0.06 


0.22 


0.16 


This work 


NDWFS Bootes 


1.6-2.9 


1.3 


8.5 


130 



^ Abbreviated field names correspond to the following: COSMOSl & COSMOS2 - two pointings 
within the Cosmic Evolution Survey (COSMOS) field, CDF-N - Chandra Deep Field North, CDF-S 
- Chandra Deep Field South, B3J2330+3927 - a high redshift radio galaxy field, GOODS-N - Great 
Observatories Origins Deep Survey North field, J2143-4423 - a galaxy cluster field, LH - Lockman 
Hole field, NDWFS Bootes - NOAO Deep Wide-Field Survey Bootes field, SDF - Subaru Deep Field 
SFLS - Spitzer First Look Survey field, SSA22 - Small Selected Area 22 field originally from lLillv et al.l 
{l991|), SXDS - Subaru/XMM-Newton Deep Survey, XMM-LSS - XMM-Newton Large Scale Structure 

Survey field. 

Exact survey area was not given in lSmith Sz Jarvi i (l2Q07h . 



Table 3 

Known Lyo; Nebulae in the NDWFS Bootes Field 



Source 


Survey^ 


Redshift Bw 


Isophotal Area^ Bw 
(arcsec^) 


Surface Brightness^ 
(mag arcsec"^) 


Notes^ 


LABd05 


Dey et al. (2005) 


2.7 


54.4 




27.0 


Recovered 


PI 


Prescott et al. - Subaru 


2.7 


27.6 




27.3 


Too Faint 


P2 


Prescott et al. - Subaru 


2.6 


71.0 




26.9 


Core-Halo 


PS 


Prescott et al. - Mayall 


1.9 


25.4 




26.9 


Core-Halo 


Yl 


Yang et al. (2009) 


2.3 


20.4 




26.7 


Core-Halo 


Y2 


Yang et al. (2009) 


2.3 


17.0 




26.8 


Core-Halo 


Y3 


Yang et al. (2009) 


2.3 


36.6 




25.9 


Core-Halo 


Y4 


Yang et al. (2009) 


2.3 


37.5 




26.6 


Core-Halo 


^ Spectroscopically confirmed Lya nebulae selected from a Subaru intermediate-band survey (Prescott et al. 2012, in 
preparation) and a Mayall intermediate-band survey (Prescott et al. 2012, in preparation) are listed along with Lya 
nebulae from the narrow-band survey by Yang et al. (2009). 

^ Isophotal area determined using SourceExtractor and a detection threshold set to the median la surface brightness 
threshold of the NDWFS Bw imaging (Section [OJ. 

Mean surface brightness within the isophotal area determined using SourceExtractor fSection|4.2|). 
^ "Too Faint" indicates that the diffuse emission was too faint to be selected by our survey. "Core-Halo" describes 
sources with a central compact continuum source surrounded by a Lyo; halo; these sources are partially or completely 
removed from the images during pipeline steps designed to remove bright stars and compact galaxies and are not 



selected by our survey. 
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Figure 1. Broad-band Bw bandpass and a generic narrow-band filter bandpass (NB; FWHM = lOOA) overplotted on a spectrum of the 
night sky. Our systematic broad-band search technique, designed to find luminous Lyo; nebulae at 2 < 2; < 3, is able to efficiently survey 
enormous comoving volumes by relying on deep, wide-field, broad-band imaging and the very dark sky in the blue. 
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Background & 
Halo Subtraction 



Bright 
Object 
Removal 



Faint 
Compact 

Object 
Removal 




Identify 
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Extended 

Objects 



<1) Source 
Extractor 



2) Source^ 
Extractor 



High Threshold 



4) Source^ 
Extractor 



Medium Threshold 



6) Unsharp^ 
V Masking 





^3) Source 
Extractor 



Low Threshold 



^ 5) Source 
Extractor 



Low Threshold 



^ 7) Source 
Extractor 

Remove 
Detections 



8) Remove 
Regions 
Around 

Bright Stars 



9) Wavelet ^ 
Decomposition 



Final Catalog 
Cleaning & 
Candidate 

Prioritization 




^10) SourceX 
Extractor / 



Select High 
Power Peaks 



11) Object^ 
Flagging 



^ 12) By-eye ^ 
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Image 



13) Size-Color 
Prioritization 



Figure 2. Flowchart outlining the individual steps of our systematic broad-band search for luminous Lyo; nebulae. 
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Figure 3. Individual steps of the broad-band Lya nebula search algorithm illustrated using the region around LABdOS (Section 12. 2^ . 
The diffuse emission of LABdOS is clearly visible (arrows) in the final frame (Step 8). 
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Figure 4. Wavelet decomposition is used to separate out sources of different size scales. The upper left panel shows the final result of 
Steps 1-8 of the search algorithm, as in Figure O The other three panels show the wavelet power maps for different size scales. The diffuse 
emission of LABdOS generates a large wavelet power peak in the final panel and is easily sele cted a s a hya nebula candidate (arrows). The 
source in the upper righthand corner of the final panel was below the 4cr threshold (Section 12. 2. 5|) and was not selected. 



A Successful Broad-band Survey for Giant Lya Nebulae I 13 



u 

CD 

u 



CD 

bo 



CD 
> 

o 



CJ 
CD 

U 



CD 
00 



CD 

> 



□ 



□□□ 



100 - 



100 - 



1^ 





-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 



Figure 5. Wavelet Size vs. Bw — R color for Lya nebula candidates selected using the morphological search pipeline. The wave let siz e 
corresponds to the size of the source in the wavelet power map, but does not indicate the true nebular size of the object (Section I2.2.7|) . 
Top: Lya nebula candidates are shown (black circles) along with artifacts from the residual halos of bright galaxies (galaxy; cyan squares) 
and detections of tidal tails and spiral arms (tidal/ ARM; green triangles). Bottom: Lya nebula candidates are shown as in the top panel 
(black circles), with first and second priority selection regions indicated (blue dashed lines). The first (second) priority region contains 39 
(40) Lya nebula candidates. LABd05 is indicated with a double star (Dev et al. 2005), and the other spectroscopically-confirmed Lya 
nebulae are shown as single stars (Prescott et al. 2009, Paper II). Artifacts from the upper panel are indicated with small black dots. The 
large filled black circle with an error bar represents the typical color of low surface brightness galaxies (LSBs; Haberzet tl et al.i i2007), and 
the histogram (plotted on a linear scale) represents the distribution of Bw — R colors for field galaxies in NDWFS, demonstrating that the 
colors of our final Lya nebula candidates are substantially bluer than typical LSBs and field galaxies. 
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Figure 6. Expected B\y — R colors for Lya nebulae at redshifts 0.1 < z < 3.1, where the colors are divided into 0.1 mag bins; the 
By[/ — R color cuts used in our survey are shown as black lines at B\y — R < 0.45 and at B\y — R < 0.65. Lya nebulae are modeled as 
flat spectrum sources (fx oc A~^) with Lya emission lines of varying equivalent widths, where the continuum at the position of Lya is 
measured using a 20A bin redward of the line. The B\y bandpass contains Lya at 1.9 ^ z < 2.9, leading to blue colors. However, due 
to the wavelength range of the spectroscopic follow-up observations, Lya nebulae can be detected as spatially extended sources even at 
lower redshifts (z > 1.6) if the continuum has the same extended morphology as the Lya, even though Lya is not contained with the 
Bw band. The canonical upper limit for Lya equivalent widths arising from stellar processes is shown (240A , dashed line: ICharlot FalU 
11993) along with the positions in redshift-EWrest space of the conflrmed Lya nebulae (Dev et al. 2005; Presc ott et al.|[20Q^ . Paper II). The 
B\Y-R colors of actual Lya nebulae can and do deviate from the colors predicted using this simple model most likely due to, e.g., differing 
continuum shapes, the composite nature of actual Lya nebula systems, the presence of embedded sources with different colors. 
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Figure 7. Sky distribution of priority 1 (filled circles) and priority 2 (crossed open circles) Lyo; nebula candidates within the NDWFS 
Bootes field. The individual NDWFS tiles are shown as blue (Bw imaging) and red (R-hand imaging) dashed rectangles. Confirmed Lyo; 
nebulae are shown as stars (Paper II), including LABdOS (double star: .Dev et al.ii2005i '). 
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Figure 8. Cumulative expected number of Lyo; nebulae assuming a constant (i.e., redshift independent) comoving volume density of 
0.5 X 10~^ — 3.6 X 10~^ Mpc~^ (the range estimated for z ^ 2.3 bv lYang et al. 2009) as a function of redshift and a 100% detection rate 
(solid lines). The dashed lines show the cumulative expected number of Lyo; nebulae under the assumption that the comoving volume 
density increases linearly to higher redshift and matches the measured range at z = 2.3. Our survey is designed to select Lyo; nebulae at 
redshifts z ~ 1.6 — 2.9 (dark grey shading); however, Lyo; nebulae at z ~ 1.2 — 1.6 (light grey shading) with diffuse blue continuum emission 
that are selected by our morphological survey will appear as continuum-only sources in follow-up spectroscopy. 




10 100 
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Figure 9. Average surface brightness versus isophotal area (above the median surface brightness limit of the NDWFS survey; 28.9 mag 
arcsec^, Icr, 1.1^^ diameter aperture) for both first (filled circles) and second (open circles) priority candidates. Confirmed Lya nebulae 
from this survey are indicated with stars (Paper II), including LABdOS (double star; Dev et al. 2005). Spectroscopically confirmed Lyce 
sources in the Bootes field from other surveys that were not selected by our broad-band search (Table [3]), primarily due to the presence 
of bright central sources, are shown as diamonds (Yang et al. 2009) and squares (Prescott et al. 2012e, in prep.). The median 5cr surface 
brightness limit of the survey is shown as the black line, with the standard deviation of all 27 pointings shown as a shaded band. 
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1) NDWFS J143004.7+353509 Priority = 2 




3) NDWFS J142846.2+330819 Priority = 3 




4) NDWFS J143436.4+330406 Priority = 2 




6) NDWFS J143153.1+351436 Priority = 2 




7) NDWFS J143055.2+352855 Priority = 1 




Figure 10. Postage stamp images of the Lya nebula candidates in the Bw, and / bands, respectively. Images are displayed using a 
log scaling within 30^^x30^^ boxes, and are labeled with the candidate number, name, and priority listed in Table]^ 



APPENDIX 
POSTAGE STAMPS 

Postage stamp images of the 79 Lya nebula candidates from our broad-band survey of the NDWFS Bootes field. 
The Bw and i?-band images shown were used for morphological and color selection during the survey; the /-band 
images are shown for reference. 
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Figure 11. Postage stamp images of the Lya nebula candidates in the Bw, and / bands, respectively. Images are displayed using a 
log scaling within 30^^x30^^ boxes, and are labeled with the candidate number, name, and priority listed in Table]^ 
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21)NDWFS J143230.2+334105 Priority = 2 26) NDWFS J142622.9+351422 Priority = 1 




22) NDWFS J142819.8+344657 Priority = 2 27) NDWFS J143055.7+340502 Priority = 2 




Figure 12. Postage stamp images of the Lya nebula candidates in the Bw, and / bands, respectively. Images are displayed using a 
log scaling within 30^^x30^^ boxes, and are labeled with the candidate number, name, and priority listed in Table]^ 
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Figure 13. Postage stamp images of the Lya nebula candidates in the Bw, and / bands, respectively. Images are displayed using a 
log scaling within 30^^x30^^ boxes, and are labeled with the candidate number, name, and priority listed in Table]^ 
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41) NDWFS J142758.6+354428 Priority = 2 




42) NDWFS J142634. 1+334954 Priority = 2 




43) NDWFS J143435.9+350710 Priority = 2 




44) NDWFS J142927.8+345906 Priority = 3 




45) NDWFS J142745.2+352949 Priority = 1 



46) NDWFS J142826.0+352116 Priority = 2 



47) NDWFS J142535.2+324934 Priority = 2 





49) NDWFS J142821.3+352001 Priority = 1 



mm 



50) NDWFS J143401.0+351229 Priority = 1 



Figure 14. Postage stamp images of the Lya nebula candidates in the Bw, and / bands, respectively. Images are displayed using a 
log scaling within 30^^x30^^ boxes, and are labeled with the candidate number, name, and priority listed in Table]^ 
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51)NDWFS J143348.9+354157 Priority = 1 56) NDWFS J143141.3+352110 Priority = 2 




54) NDWFS J142522.3+325424 Priority = 2 59) NDWFS J143412.7+332939 Priority = 1 




Figure 15. Postage stamp images of the Lya nebula candidates in the Bw, and / bands, respectively. Images are displayed using a 
log scaling within 30^^x30^^ boxes, and are labeled with the candidate number, name, and priority listed in Table]^ 
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61) NDWFS J143438.7+350420 Priority = 2 
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62) NDWFS J142832.7+354105 Priority = 1 



63) NDWFS J143335.8+344503 Priority = 1 



[■-■■ ■ • ]■ ' • 1 
1 • • 

i; : *: ^: i: i: : . 




64) NDWFS J143203.8+351855 Priority = 2 







65) NDWFS J143207.2+343101 Priority = 3 

|: : ■ : ■': I : 




66) NDWFS J142539.9+344959 Priority = 1 



67) NDWFS J142819.7+325449 Priority = 1 



I. 



68) NDWFS J142443.9+344834 Priority = 1 



69) NDWFS J143202.6+351904 Priority = 1 




70) NDWFS J142753.8+341204 Priority = 1 




Figure 16. Postage stamp images of the Lya nebula candidates in the Bw, and / bands, respectively. Images are displayed using a 
log scaling within 30^^x30^^ boxes, and are labeled with the candidate number, name, and priority listed in Table]^ 
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Figure 17. Postage stamp images of the Lya nebula candidates in the Bw, and / bands, respectively. Images are displayed using a 
log scaling within 30^^x30^^ boxes, and are labeled with the candidate number, name, and priority listed in Table]^ 
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81) NDWFS J142735.5+342332 Priority = 1 




82) NDWFS J142449.8+324743 Priority = 1 




83) NDWFS J142732.5+341213 Priority = 1 




Figure 18. Postage stamp images of the Lya nebula candidates in the Bw, and / bands, respectively. Images are displayed using a 
log scaling within 30^^x30^^ boxes, and are labeled with the candidate number, name, and priority listed in Table [T] 



